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Abstract. The correction to the string junction three-quark potential in a baryon due to the proper moment 
of inertia of the QCD string is calculated. The magnitudes of the string corrections in P wave heavy baryons 
are estimated. 

PACS. 12.38.-t Quantum chromodynamics - 12.40.Yx Hadron mass models and calculations 



1 Introduction 



The formation of an extended object, the QCD string, 
between the color constituents inside hadrons plays cru- 
cial role in understanding their properties [T] . In a quark- 
antiquark system, the string potential V(r) = a r, where 
cr is the string tension, represents only the leading term 
in the expansion of the QCD string Hamiltonian in terms 
of angular velocities. The leading correction in this expan- 
sion is known as a string correction. This correction works 
well at small orbital excitations, when there are no long 
strings. So far, the string correction was calculeted for the 
orbitally excited heavy-light mesons [2] and hybrid char- 
monium states [3] . The opposite case of the large angular 
excitations is way more advanced. 

In a baryon, the leading term of the string potential is 



W(ri, r 2 , r 3 ) = ar, n 



(1) 



where r m i n is the minimal string length corresponding to 
the Y-shaped configuration. In this paper, we generalize 
the result of Ref. [2] and derive the string corrections to 
the string junction potential ([T]). As an application, we 
consider the P-wave excitations of heavy baryons. 

The plan of the paper is as follow. In section [21 we 
briefly review the derivation of the string potential as ap- 
plied for low orbital excitations of quark-antiquark sys- 
tems. Section 3 is devoted to the generalization of this 
derivation for the baryon Wilson loop. Section 4 presents 
the numerical results for heavy baryons. Finally, Sect. 5 
summarizes our conclusions. Some details of the calcula- 
tion of the baryon string correction are explained in Ap- 
pendix. 
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2 String QQ potential 

First, we recall how the string interaction arizes in a quark- 
antiquark system. We create a gauge invariant quark- 
antiquark state at time T — which is annihilated at a 
later time T. It is usually assumed that the string poten- 
tial appears from the minimal area law asymptotic for an 
isolated Wilson loop 



< W > = 
= ^TrPexp 



ig f dz n A n 



exp(-CTS , min ), (2) 



where the left-hand side is the average over the stochastic 
QCD vacuum and S m i n is the area of the minimal 
surface bounded by the contour C swept by the quark 
and antiquark trajectories, 



string 



lim — In < W > 

T-s-oo T 



(3) 



In terms of functions w a (t, /3) which map a position in 
the parameter space (r, ft ) to a point in space-time the 
minimal area S^in in Eq. @ is written as 



oIt / d(3^(ww>y 



(4) 



where = dw u /dr, w'^ = dw u /dj3. We use the instanta- 
neous approximation t\ = t% = t and the laboratory gauge 
t = t. Choosing for w^(t, 0) the linear form corresponding 
to straight-line string 



w^t, (3) = p x ln (t) + (1 - (3) x 2n {t), 
one obtains for the Lagrangian 



(5) 



L = 



nil 



™2_ 

2^2 



Mi + M2 



• / dPy/1 - [n 
Jo 



x 08fi + (1 - /8)f a )] 



(6) 
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where r = r\—r2, n = r/r, j3 is the remaining internal pa- 
rameter along the string, the dots denote derivatives with 
respect to the physical time t, and are the bare quark 
masses. The angular velocities in ^ describes the con- 
tribution of the proper inertia of the rotating string. The 
Lagrangian © is written in the einbein field form [5] ; the 
einbeins fi%, fi2 are introduced to simplify the treatment of 
relativistic kinematics. In the center-of-mass (CM) system 
r r 
" 2' 



ri = 2 



r 2 



Pi ■ P2 

ri = — , r 2 = — , p 1 
Mi Mi 



P 2 = P 



(7) 



The general procedure for the transition from the La- 
grangian ^ to the quantum Hamiltonian suggested in [JJ 
is rather complicated due to the presence of square roots 
in However, if one is interested in the low- lying part 
of the spectrum, the potential-type regime can be consid- 
ered. The leading term corresponds to the I independent 
term, in the expansion of the string potential in powers 
of I 2 , where I — r x (— iV r ) is the orbital momentum op- 
erator. The first correction, of order I , is known as the 
string correction. 

The Hamiltonian function corresponding to © reads 



P A 



2/1; 



V f 



string 



(r), 



(8) 



where 
V, 



string 



(r) 



nx Ml-(l-ffl£ 
Mi M2 



(9) 



Expanding the square root and making the substitution 
p — > — i V r one obtains an approximate expression for the 
string potential [2] 



V. 



strin 



err 



or 




(1 



M2 



mI 



1-/3 

M2 



M1M2) I 2 



d;5 



err — - - - — (10) 

The last term in Eq. (ITU)) is the string correction. Its sign 
is negative so that the contribution of the string correc- 
tions lowers the energy of the system, thus giving a nega- 
tive contribution to the masses of orbitally excited states, 
leaving those with I = intact. 



3 The string baryon potential 

Now let us apply a similar consideration to the baryon 
states. In tricolor QCD, three quarks are produced at 
the origin and then fly, eventually creating strings of the 



4^ 



a'b'c' 




Fig. 1. The baryonic Wilson loop. The quarks are located at 
positions r\,r2 and r?,. 



chromo electric field. If the three quarks are at positions 
7*1, r2, and r3, the strings meet at an interior point, known 
as the string junction point, where the angles between the 
flux tubes are 120°, provided that none of the angles of 
the triangle formed by the three quarks exceeds 120°. Oth- 
erwise, the configuration of minimal total length is made 
of two flux tubes meeting at the third quark location. In 
the general case time evolution produces a three-bladed 
Wilson loop shown in Fig. [1] 

For a baryon with the genuine string junction point the 
calculation of the string correction is a very cumbersome 
problem. The calculations are greatly simplified, however, 
if the string junction point is chosen to coincide with the 
CM coordinate R cm . In this case, the complicated string 
junction potential is approximated by a sum of the one- 
body confining potentials. The accuracy of this approxi- 
mation for the unperturbed three-quark Hamiltonian was 
tested in [6]. In particular, for the P-wave baryon states 
the accuracy is better than 1% . 

Letting R cm = 0, we arrive at the following general- 
ization of Eq. (|10p for a baryon 



V. 



CM 



string 



E 



d(3 t 



•l- 



E 
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2 2 



I 



1 - - 



6 nf rf 



(11) 



where li = r.i x p i are the orbital momentum opera- 
tors of the individual quarks. Each term in Eq. (|TT|) is 
associated with the string beginning at the center-of-mass 
coordinate and ending at the quark position r^. As a re- 
sult, for the string correction we get 



AM. 



string 



(12) 
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where \P is an eigenfunction of the unperturbed Hamilto- 
nian. 



4 The string correction 

The calculation of the matrix element (TT21 is more con- 
veniently performed in terms of the Jacoby variables p, A 
rather than quark coordinates r,. The baryon wave func- 
tion depends on the three-body Jacobi coordinates 




\.. = JBhl ( ^ iTi + till _ ,., 



(13) 



Mo V Mi + Mj 
where fiij and are the appropriate reduced masses: 



M/ 



MiMj 



Mi 



Mj 



ftij, k 



M 



(14) 



ilf = /ii + fi2 +M3j an d Mo is an arbitrary parameter with 
the dimension of mass, which drops out in the final ex- 
pressions. The hyperradius R 2 = po {pfj + \j) does not 
depend on po and is invariant under permutations. There 
are three equivalent ways of introducing the coordinates 
(fT3|) . which are related to each other by linear transfor- 
mations with the Jacobian equal to unity. 

The quark orbital momenta li in Eq. (Ill[) can be sim- 
ply constructed from the orbital momenta corresponding 
to different Jacoby variables. It is convenient to choose 
for each term containing I 2 in (|lip an appropriate set of 
Jacoby coordinates pj k , then we obtain 



M?^ 



where ijk = 123, 231, 312, and 




Mo MiA ifc 



H) 



8 



d\jk 



(15) 



(16) 



In what follows we consider the nnQ baryons, where 
n stands for the light quarks u or d, and Q stands for the 
heavy quarks c or b. We use the basis in which the heavy 
Q-quark is singled out as quark 3 but in which the nn 
quarks are still antisymmetrized. The similar basis is used 
for the nsQ baryons. These basis states diagonalize the 
confinement problem with eigenfunctions that correspond 
to separate excitations of the light and heavy quarks. We 
refer to these excitations as the p- and A-excitations, re- 
spectively. The nonsymmetrized nnQ and nsQ bases usu- 
ally provide a much simplified picture of the states. The 
physical P-wave states are neither pure SU(3) states nor 
pure p or A excitations but linear combinations of all states 
with a given J. Most physical states are, however, closer 
to pure p or A states than to pure SU(3) states [7]. 

For definiteness we now discuss the simplest case of 
P-wave states. The angular functions of these states can 



be of two types. In Cartesian coordinates they are given 
by 



Y, 



6 Pig 

7T 3 R 



«31 ~FT + P31 



R 



R 



P23 
R 



a23-Tf + $23 



A23 

~R 



(17) 



and 




331- 



P31 
R 



CH31- 



A31 
~R 



a P23 , „ ^23 
P23-^- + "23 



(18) 



where the mass dependent coefficients aij and /3V,- are de- 
fined in (|A.4I) . (|A.5I) . The functions Y p> \ are the eigen- 
values of the six-dimensional angular operator 

Z 2 Y p , A (#,n p ,n A ) = - K(K + 4)Y PiX (e,n p ,n x ),(19) 

where 



o 2 



8 l 2 
4cot6> — £- 



00 



cos z 



(20) 



and K = 1 is the grand orbital momentum. The functions 
\P V in (IT2|) are written as 



u„{R) 



Y„(0, n p ,n x ), v = p, A, 



(21) 



^5/2 ' *v>^' 

where we have introduced the reduced radial functions 
u v (R). The further calculations are elementary and lead 
to the final resuliQ 

64 a 



1 



45tt VM (pi + p 2 ) 



VM2 + M3 

3/2 
, Mi' 



M2 



AM 



64a 



VMi + M 3 
3/2 
M2 



M3 



"Mi 7p . 



(22) 



457T M 3 / 2 (pi + p 2 ) 

/'i +M2I 5 , /M 2 + M3 
Mi 



where 



M3 



U v {x) 



Mi + M3 



M2 



c?x, v — p, X, 



7a, (23) 



(24) 



x = y/JiQ R, and the functions are normalized as 

u 2 (x)da; = 1 (25) 

Eqs. (|2"2"]) . (|2"3")) define the magnitude of the string correc- 
tion for the P-wave baryon states. The generalization to 
the D-wave excitations is straightforward. 



1 Some details of the calculations are presented in Appendix 
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5 Numerical results 

Now we are in position to estimate the magnitudes of 
string corrections. Here we concentrate on heavy baryons. 

Our calculations of the unperturbed baryon masses 
and wave functions are performed using the effective Hamil- 
tonian derived in the Field Correlator Method [8]. For 
the different technical aspects of these calculation see refs. 
[H]-[TT]. As an unperturbed Hamiltonian we consider the 
Hamiltonian specified in Ref. [B] 

3 2 

^^iMr.r.r^l/c, (26) 

i=l ^ 

where the long-range confining force (JT|) is augmented by 
the short-range Coulomb potential, 

As in Eq. (j6j , an approximate einbein field method is used 
for derivation of Eq. ([2"6"]l : einbeins fa are treated as c- 
number variational parameters. The eigenvalues Eo(mi, fa) 
of the Hamiltonian (|26|) are found as functions of the bare 
quark masses rrii and einbeins fa, and are finally mini- 
mized with respect to the fa. With such simplifying as- 
sumptions the spinless Hamiltonian Ho takes an appar- 
ently nonrelativistic form, with einbein fields playing the 
role of the constituent masses of the quarks. 
The zero-order baryon mass Mq is 

M ° = E(|f +y) +Eo(fa)+C (28) 

In Eq. (|28p a constant C is the sum of the nonperturba- 
tive self-energies of the quarks propagating through the 
vacuum background field which provides an overall nega- 
tive shift of Mo required by phenomenology [12.. Finally, 
Eo(jj,{) is an eigenvalue of the Schrodinger equation 

d 2 u v (R) 
dB 2 + 

2 [E - ^ - V$(R) - V£(R)\ u v (R) = 0, (29) 

where Vy(R), Vq(R) are the averages of the string junc- 
tion and Coulomb potentials, respectively, over the six- 
dimensional sphere £l§ . Explicit expressions of these quan- 
tities can be found in Ref. [TU]. We use the same param- 
eters, a = 0.15 GeV, a s = 0.39, m„ = 7 MeV, m s = 185 
MeV, m c = 1.36 GeV, and m b = 4.71 GeV, that were 
used in Ref. fT3] for the analysis of the ground states of 
heavy baryons. Thus the model is totally fixed. 

The wave function factors j u in Eqs. ([22)1 . ([23")) which 
define the magnitude of the string correction do not es- 
sentially depend on the baryon flavor nor on the type of 
excitation. This property mentioned for the light baryons 



in Ref. [TT] is also valid for heavy baryons. E.g. for all con- 
sidered baryons we get j p , 7a ~ 0.30-0.31 GeV 1 ' 2 , which 
are practically the same values as were found for the light 
baryons. The relative size of the string corrections for dif- 
ferent baryons is mainly defined by the constituent quark 
masses fa that are themselves the calculable quantities. 

Note that our technique does not take into account 
chiral degrees of freedom, which are essential for light di- 
quarks in the Aq baryons, so we consider the Eq and Eq 
where the chiral effects are expected to play the minor 
role. The result of the calculation of the P-wave states is 
given in Table [T] In this table we present the dynamical 
quark masses /x n , fa and fiQ (columns 3-5) for various 
heavy baryons. The latter are computed solely in terms of 
the bare quark masses, a and a s and marginally depend 
on a baryon. We also display the eigenvalues Eo of the 
Schrodinger equation (column 6), the zero-order masses 
calculated without the string corrections from Eq. ([28)) 
(column 7), the string correction itself (column 8) and the 
total baryon masses (column 9). Note that the A excita- 
tions are typically ~ 100 MeV heavier than the p excita- 
tions. 

Because in this paper we neglect the spin interactions, 
it would be premature to compare literally the results 
of our piloting calculations with experimental data for 
the Eq and Eq excitations. Presently, the Particle Data 

Group [14] lists only two | c baryons, 2^(2800) and 
£■(,((2790). Spin and parity of these states have not been 
measured but were assigned on the basis of quark model 
predictions. The P-wave b baryons have not been detected 
yet. In Table 2 we compare our predictions for the A exci- 
tations with the results of several quark models [IS] - [IS] . 
It is seen that our results match the region of the predic- 
tions well. 



6 Conclusions 

In this paper, we have shown how the dynamics of the 
strings with quarks at the ends influences the masses of 
the baryon orbital excited states. Note that the string cor- 
rection is totally missing in quark models based on the 
relativistic equations with local potentials. We extend our 
approach to baryons containing charm and bottom quarks 
and calculate the string correction to the P-wave states of 
the nnQ and nsQ baryons. Our main result is given by 
Eqs. ([22)) . ([23")l . The numerical results are given in Table [1] 
For each baryon, we have calculated the dynamical quark 
masses /i,-, the baryon masses with the self-energy correc- 
tions in Eq. (125)) and the string corrections. The latter 
provide an extra piece of the effective interquark poten- 
tial, which is entirely due to the string-type interaction in 
QCD. We emphasize that no fitting parameters were used 
in baryon spectrum calculations. The estimated string cor- 
rections give a negative contribution to the baryon masses 
of the order 30-50 MeV. Our comparative study provides 
deeper insight into the quark model results for which the 
string interactions encode the underlying QCD dynamics. 
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Appendix 

In this appendix, we present for an illustrative purpose 
the explicit calculation of the first term with i = 1 in Eq. 
(fTTj) . Consider as an example the p excitation. Eq. (fT5|) 
implies that 



/ 2 



mfi 



\%) = 



M2 + M3 




1 

A23 



(A.l) 



M y y Mo Ml 
At this point it is convenient to restore the quark numera- 
tion in the Jacobi coordinates. We use the transformation 
law 

«31 #51 \ ( P31 

■/?31 "31/ \^31. 



with 



«3i = 



^31 



"23 



2.'! 



«23 #23 \ ( P23 
^ /^23 "23 / \ A 2 3 

M2M3 

(Ml +M2)Ol +M3) 

(Mi + M2XM1 + M3) ' 

/ MlM3 

V (Ml + M2)(M2 +M3) 

/ m 2 M 

V (Mi + M2XM2 + M3) 



(A.2) 
(A.3) 



(A.4) 



(A.5) 



Substituting 3^(12; 3) = a 23 ^ p (23, 1) + /3 23 <fp(23, 1) and 
noting that 



P23 
; R 



0, 



-^23 

~R 



~R 



(A.6) 



we obtain for the matrix clement in the r.h.s. of (IA.1 



1 



\ 2:j 



\*p) = 



= 2f3 23 I ("aafea + Pa 3 A a3 )A a 3 df? 



i? 2 



di? = 



A 2 3 
647T 2 i? 



45 



023 



i? 



di?, 



(A.7) 



which corresponds to the first term in Eq. (12"2"1) . An anal- 
ogous calculation with the obvious changes leads to the 
second term in (|2"2"|) . The string correction to the A exci- 
tation can be evaluated in a similar way. 
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Table 1. Heavy baryons with L = 1 



Baryon L Q fi n ^ s Eo M AM stv i ng M B 
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4764 
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Table 2. Theoretical predictions for masses of heavy baryons in various quark models 
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